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Abstract 

We study nonperturbative corrections to the inclusive rare decay B 
Xsi~^i~ by performing an operator product expansion (OPE) to 0{l/m'^) 
The values of the matrix elements entering at this order are unknown and 
introduce uncertainties into physical quantities. We study uncertainties in- 
troduced into the partially integrated rate, moments of the hadronic spectrum, 
as well as the forward-backward asymmetry. We find that for large dilepton 
invariant mass > M^, these uncertainties are large. We also assess the 
possibility of extracting the HQET parameters Ai and A using data from this 
process. 
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I. INTRODUCTION 



Rare B decays mediated via flavour changing neutral currents have received much at- 
tention because of their sensitivity to physics beyond the standard model In the standard 
model these decays occur via penguin and box diagrams with virtual electroweak bosons 
and up-type quarks in the loops. Because of the large top quark mass, the contribution with 
a top quark in the loop dominates. At energy scales below the mass of the top quark and 
W boson, it is convenient to switch to an effective theory where the top quark and the weak 
bosons have been integrated out of the theory. The 6 — > s transition is then mediated by 
the effective Hamiltonian [1] 

AC ^° 

neff = -^v,,vzj:Q{^i)o,(^i), (1) 

where the operators are commonly defined by 

01 = {sLa1l^bLa){cL0l^CL(3), 

02 = {sLaliJ>LI3){cLI3YcLa), 

03 = {sLal^LbLa) (QL^^^QLo) , 

q=u,d,s,c,b 

q=u,d,s,c,b 
q=u.,d,s,c,b 
q=u,d,s,c,b 

O7 = -:^^SaCr^^{mbR + msL)baF'"', 
e2 _ 

Olo = r^s^rLbJ^^^.i. (2) 

Here L/R — ^(l± 7^) are the usual left and right handed chiral projection operators. The 
values of the Wilson coefficients Ci(mb) have been calculated in the next to leading log 
approximation [2,3] in the standard model and are given in Table I. 



Ci 


C2 


C3 


C4 


C5 


Ce 


C7 


Cg 


Cio 


-0.240 


1.103 


0.011 


-0.025 


0.007 


-0.030 


-0.311 


4.153 


-4.546 



TABLE 1. The Wilson coefficients Ci{mh) in the next-to- leading log approximation. 



Physics beyond the standard model will generally introduce new contributions to the 
loop and will therefore modify the values of these coefficients [4]. Measurements of the 
Wilson coefficients may therefore indirectly constrain new physics scenarios. For example. 
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the decay 6 ^ S7 is proportional to jCyp, and the recent measurements of the branching 
ratio B — > K*j [5] and inclusive rate B Xg'-f [6] have placed constraints on models of 
physics beyond the standard model which modify the magnitude of C7 [7]. 

The decay b — > si'^i" is suppressed, relative to 6 — > S7, by an additional factor of 
the electromagnetic coupling constant and has not yet been observed [8]. It has, however, 
the appeal of being sensitive to the signs and magnitudes of C7, Cg, and Ciq, making it 
a potentially more powerful probe than 6 ^ 57 of beyond the standard model physics. 
Experimental studies of this process impose cuts on the available phase space. This is 
primarily due to the necessity of removing the resonance from B — > [J /ip,ip')Xs with the 
( J/V', ■00 decaying into two leptons. We incorporate representative cuts into the theoretical 
analysis. 

Using an operator product expansion (OPE) several observables of the inclusive decay 
B —>■ Xst^l~ have been calculated including the leading non-perturbative corrections [9,10]. 
In this framework these leading corrections arise as matrix elements of dimension five oper- 
ators, suppressed by two powers of the h quark mass, and are conventionally parameterized 
by two quantities, Ai and A2. A third parameter A enters through the difference of the h 
quark mass and B meson mass 

m, = Ms-K + ^^^P^ + --: (3) 

Whereas A2 can be determined from the B* — B mass sphtting, A2 = (M^. — M^)/4 ~ 
0.12 GeV^, no such simple relation exists for Ai. It has been estimated using various methods 
to lie in the range 0.1 GeV^ < (-Ai) < 0.6 GeV^ [12]. 

In a previous paper we extended the analysis of the total rate for B XsC.^C~ to one 
order higher in the OPE [13]. The dimension six operators arising at this order can be 
parametrized by six quantities, commonly labelled pi_2 and 7^_4, all of which are unknown. 
We found that the uncertainties introduced by these six parameters can be significant, 
depending primarily on the actual values of the matrix elements and the amount of accessible 
phase space. In this paper we give the details of that analysis and also present calculations 
for the forward-backward asymmetry and moments of the hadron invariant mass spectrum 
at 0{l/ml). As in our previous analysis, we neglect perturbative effects and effects due to 
the finite mass of the s quark, which have been considered elsewhere [10]. 

It has also been proposed that, rather than use this decay to search for new physics, it 
might instead be used to extract the parameters A and Ai through a measurement of its 
hadronic invariant mass moments [10]. We estimate the uncertainties in this extraction due 
to the unknown matrix elements of the dimension six operators. 

This paper is organized as follows. In section II we briefiy introduce the formalism used 
to calculate the nonperturbative corrections, and we present the results for the decay rate in 
section 111. In section IV we calculate the forward-backward asymmetry of the lepton pair. 
We then proceed in section V to calculate moments of the hadronic invariant mass spectrum 
and estimate uncertainties in extracting A and Ai from these moments. Finally we discuss 
the results and state our conclusions. 
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II. OPERATOR PRODUCT EXPANSION AND KINEMATICS 



The procedure for calculating nonperturbative contributions to heavy hadron decays has 
been thoroughly discussed in the literature [14,15], and we present here only a brief outline 
of the technique. The differential rate is proportional to the product of a lepton tensor L^,^ 
and a hadron tensor W^" and for the process in question it may be written as 



2Mb 27r2 



(4) 



where 11 denotes the three body phase space. The spin-summed tensor L^,, for massless 
leptons is 



(5) 



The hadron tensor W^'^ is related via the optical theorem to the imaginary part of the 
forward scattering matrix element W^^ — (— I/tt) ImT'^'^ where 



-i j d^xe-"!-'^ {b |r{J^(^)'(x), J^^'^\Q)}\ b) . (6) 
In this equation denotes the current mediating this transition, and is given by 



(7) 



where q = (p+ is the dilepton momentum^. In accordance with convention, we have 

defined two effective Wilson coefficients: C^^^ = C7 — C5/3 — Cq and C^^^ . The latter 
contains the operator mixing of Oi^e into O9 as well as the one loop matrix elements of 
Oi_6,9 [2,3]. The full analytic expression for C^^^ is quite lengthy and may be found in [3]. 

Since in the decay of a 6 quark the momentum transfer to the final state parton is large, 
the time-ordered product (6) can be expanded in terms of local operators [14,15] 



i f d^xe-"'''T{j\x), J(0)} ~ — 



2mh 



Ami 



O2 + 



1 



8m] 



■O3 + . . . 



(8) 



where On represents a set of local operators of dimension d = (3 + n) , each operator con- 
taining n derivatives. For a generic current Jfj, the expressions for these operators arc quite 
lengthy. The complete set of operators for 0? < 5 [16] and d = 6 [17] appear in the literature. 
In this study we include operators up to and including dimension d — 6. 

The standard Lorentz decomposition for the forward scattering amplitude is 



-Tl Qi^^ + T2 Vf, Vy 



3 ieM-/3 f + T4 r r + v + rvn, 



(9) 



^ Notice that this current reduces to the {V — A) current when Cj-' = 0, Cg^ = 1/2, and 
Cio = —1/2. This provides some useful cross-checks with known results for semileptonic B decays 
[18]. 
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where v'^ is the four-velocity of the initial b quark ph = rribV. Since in this paper wc treat 
the final state leptons as massless {£ — e,ii), the form factors T4_5 do not contribute to 
observables. 

It is clear from (6) and (8) that to calculate these form factors we must take matrix 
elements of the operators On- Matrix elements of dimension four operators vanish at leading 
order in the l/rub expansion [14] and matrix elements of dimension five operators may be 
parameterized by Ai and A2 [19] 

{B{v)\KUD^tD,K\B{v)) = MbTt |rP+ Qai(5^, - v^v,) + ^As^a^.) P+} , (10) 

where P+ = |(1 + and F is an arbitrary Dirac structure. 

Finally, the dimension six operators may be parameterized by the matrix elements of 
two local operators [18,20] 



{B{v)\hyiDJD^,iDfih^\B{v)) = ipi {gafs -VaVf})v^, 



2M, 



2M, 

{B{v)\hyiDaiDu,iDp-fs-f5K\B{v)) = -p2ievapsv''v^. 



(11) 



and by matrix elements of two time-ordered products 



{B{v)\K{iDfhyi / d^x / dt Ci{x)\B{v)) + h.c. 



2Mf 



-l-{B{v)\K^{-ia,,)G^^Ki J d'x fdtCi{x)\B{v)) + h.c. 



rrib 
% + 3% 
nib 



(12) 



arising from a mismatch between the states \B{v)) of the effective theory and \B) of the 
full theory. The contributions from Tl_4 can most easily be incorporated by making the 
replacements [18] 



Ai ^ Ai + 



Ao 



A2 + 



rrib 
SrUh 



(13) 



in the parton level results. In addition, as we will show later there is a contribution to the 
total rate from the dimension six four-fermion operator 



^bs 



167r^ 



the matrix element of which we define as 

^ -{B\0 



2M, 



{V-A) 



(14) 



(15) 



The form factors up to 0{l/m^) have appeared in the literature [9]. The 0{l/ml) contribu- 
tions to the form factors proportional to pi_2 are presented in Appendix A. The dependence 
on 7^_4 is obtained by making the replacements (13) in the 0{l/m1) form factors. 
The triple differential branching ratio is given by 
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2Bn 



du dsdv ■ q 2Mb 



Im {2 s (T/^(t; ■q,s)+ T^{v ■ q, s 



+ {{v ■ qf - ^ - X j (^^''(^ • I s) + T^{v ■ q, s)) 
+ us (Ti{vq,s)-T,''{vq,s))}, 



(16) 



where we have defined kinematic variables v ■ q 



V ■ q, s 



m 



(f, and u 



1 

nii 



P-Y — {Pb — p+Y]. In terms of these leptonic variables the limits of phase space 



are given by 



\js + -<vq<^s + - 

-u{s,'ms) <u < u^s.rhs) 
Am^ < .s < (1 — rhsf , 



(17) 



where u{s, rhs) = y [s — (1 + ^s)^] — (1 — iri'sY]- 

For the calculation of the hadron invariant mass moments it will be convenient to express 
the phase space in terms of the parton energy fraction Xq = Eq/miy and the parton invariant 
mass fraction sq = p'^/ml. They are related to the leptonic variables introduced above via 



V ■ q — 1 — Xo 

S — 1 + Sq — 2xo. 



(18) 



The phase space can then be expressed as 



- 2^x1 - Sq <U< 2yjxl - 



So 



"^s < So < Xq 



(19) 



Since the form factors Tj are independent of this first integration is trivial and we arrive 
at 



d'^B 



16 -Bo 



dxo dso 2Mb \ tt 



1 



xl - Solm< 



1-2x0 + So) [Tl^ixo, so) + T^\xo, sq 



+ 



Xo - So 



{T^{xorso)+Ti{ 



Xo^ So, 



(20) 



In the above expressions we use the same conventions as in [9,10] and normalize the 
B — > Xsl'^l~ branching ratio to the semileptonic branching ratio 



dB{B ^ X,£+r ) = Bsi 



dT{B xj+r) 
r{B^XJue) ' 



(21) 



introducing the normalization constant 
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2 1 


167r2 




2 


f{rhc)K{rhc) 



^o = i3.„4V^-— -— ^. (22) 



In this expression f{mc) is the well-known phase space factor for the b — > ceu parton decay 
rate 

f{fhc) = 1 - Sfhl + Sfhl -ml- 24m^ log fhc (23) 

and K{mc) includes the 0{as) QCD corrections as well as the nonperturbative corrections 
up to 0{l/ml) 



where 

girhc) 



f{rhc) 



hi{mc) = Ai + — ^ (-9 + 2Aml - 12m\ + 12m\ - 15m^ - 72m^logm< 



hi{m^ = -h^ (77 - 88m^ + 24m^ - 8m^ - 5m^ + 96 log me + 72m^logme) 

+ — ^ (27 - 72m^ + 216m^ - 216m^ + 45m^ + 216m^logmc) (25) 
The analytic expression for the perturbative function Ao(mc) can be found in [21]. 



III. THE PARTIALLY INTEGRATED BRANCHING RATIO 



An interesting experimentally accessible quantity is the dilepton invariant mass spectrum. 
Evaluating the u integral in (16), and doing the integral over v-g by picking out the residues, 
we find for the dilepton invariant mass spectrum 



l(l-l«l + 2l)(2 + ^) + (l-15.^+10.3)(^-^) 



lOs^ + 23s^ - 9s^ + 13s + 11 pi 



!S3 



S2 



9(1 



ml 



+ 



.)^(2 + .)(2+AL)^4(-6-3. + 5.3)(- -) 
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P2 



4(5s^ + 19s3 + 9s2 - 7s + 22) pi 



9(1 - s) 



mi 



effi2 



+ 



rrib 



P2 

ml 



4(3s3 - 17s2 + S-3) pi 



3(1 -s) 



mt 



Re{cf{s))cf 



^0(1 — s) 

2> ml ^ ' 



Cio 



cf{s) + 2Cf^' 



(26) 
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The dependence on '7i_4 can be obtained by making the replacements (13) in (26). In this 
expression we have taken the hmit rhg 0. The corresponding expression with full rhs 
dependence is given in Appendix B. 
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FIG. 1. The differential decay spectrum. The solid line shows the parton model prediction, the 
long-dashed line includes the 0{l/ml) corrections and the short-dashed line contains all corrections 
up to Oil/ml). 



A plot of this distribution is shown in Fig. 1, 
nonperturbative matrix elements 



where we have used the values for the 



Ai = -0.19 GeV^ A2 = 0.12 GeV^ 



(27) 



For the matrix elements of the dimension six operators we use the generic size (Aqcd)'^ ^ 
(0.5 GeV)^ as suggested by dimensional analysis. The vacuum saturation approximation 
[22] predicts pi > 0, as shown, and we find the displayed spectrum is fairly insensitive to 
the sign of the other dimension six matrix elements. One immediately notices divergences 
at both endpoints of this spectrum. The divergence at the s ^ endpoint is due to the 
intermediate photon going on-shell and is a well known feature of the decay B — > Xs£'^£~ 
[1,9]. In this limit one expects this spectrum to reduce to the B — > rate with an on-shell 
photon in the final state, convoluted with the fragmentation function giving the probability 
for a photon to fragment into a lepton pair. This correspondence is explicitly verified by the 
analytic form of the divergent term 



1 dB 
Bo ds 



32 



€ff|2 



1 + 



Ai - 9 A2 1 Ipi - 27p2 Ti + 37^ - 3 {% + 3% 



2ml 



dmf 



+ 



2ml 



(28) 



where the term multiplying 1/s is proportional to the total rate for B — >■ [17]. As 
mentioned above, experimental cuts require us to stay away from this endpoint, and therefore 
automatically regulate this divergence. 

The divergence at the s ^ 1 endpoint is entirely due to the 1/m^ operators as can be 
seen from Fig. 1. In this case the analytic form of the divergent term is 



]_ dB 
Bq ds 



32 
3 ml 



fi2 

^10 



+ {2Cf + cf{s)f 



Pi 



1-s 



(29) 
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This leads, Tipon integration, to an unphysical logarithmic divergence in the expression for 
the total rate that is regulated by the mass of the s quark. (Of course, it is only consistent 
to include the mass of the s quark in the upper limit of integration if one uses the spectrum 
with the full dependence as given in Appendix B.) This divergence can be understood 
by considering a similar effect in the semileptonic decay B — > XciPe [23]. In that context, 
the origin of this divergence can be clarified by performing an OPE for the total, rather 
than the differential, rate [23,24]. Including dimension six operators in this OPE one finds 
a four fermion operator of the form 

— ^b-f'^LccYLb {g^„ - v^Vy) (30) 

contributing to the rate. In [23] the matrix element of this operator was calculated at leading 
order in perturbation theory by integrating out the c quark and its contribution to the total 
rate was found to be pilog(mc)- To calculate this matrix element it was essential that the 
mass of the c quark be large compared to the QCD scale Aqcd- Consequently, for the decay 
B — > Xsi'^i" where the same operator with s quarks rather than c quarks appears, these 
methods are not applicable because the s quark is too light. Including higher orders in 
perturbation theory the matrix element of the four fermion operator contains logarithms of 
the form log"'''^(ms) which are of order unity, making a perturbative calculation of this 
matrix element impossible. Thus, a seventh non-perturbative matrix element /i defined in 
(15) is required. It contributes only at the s — > 1 endpoint of the spectrum and cancels the 
logarithmic divergence proportional to pi log(m5) in the total rate 



f - f - (C!, + (2 Cf + Cf{s)r) S{l-s){p^ log(m.) - A) . (31) 



Another noticeable feature in the dilepton invariant mass spectrum is the cusp due to the 
cc threshold. Near this value of s the methods we used to calculate the physical spectrum 
fail because of long distance contributions from the resonant decay B XsJ/i/j, where the 
J/ip subsequently decays into two leptons. Experimentally one deals with this resonance 
region by simply cutting it out. Thus, to compare reliably to experiment we should include 
such a cut in our calculation. Defining the partially integrated branching ratio by 

we plot the contribution of the individual matrix elements relative to the leading order 
parton result in Fig. 2. For the generic size pi ~ (Aq^d) used in this plot, the contribution 
from pi is of the same size as the contribution from dimension five operators. This implies 
that including the 0{l/ml) corrections for this decay does not significantly decrease the 
nonperturbative uncertainties. We see that the nonperturbative contributions become more 
dominant as the accessible phase space is decreased^. For % ~ .75 the uncertainty from 



^We emphasize that the sizes of the pi contributions shown here should not be taken as accurate 
indications of the actual size of the corrections, but rather as estimates of the uncertainty in the 
prediction. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 

1 



FIG. 2. The fractional contributions to By^ with respect to the parton model result from the 

0{l/mD operators. The solid, dashed and dotted lines correspond to the contributions from A2, 
pi and p2, respectively. The contribution from Ai is too small to be seen. The two vertical lines 
illustrate the positions of the J/il) and the ij)' resonance. 



the pi matrix element is of the same size as the parton model prediction. This is a clear 
signal that the OPE is no longer valid if the phase space is restricted to be too close to the 
endpoint .s = 1. This breakdown of the OPE close to the cndpoint is a well-known feature 
encountered in this approach to the study of inclusive decays [26]. Unfortunately, in this 
endpoint region a shape function does not exist and an alternate approach, such as heavy 
hadron chiral perturbation theory for exclusive final states, must be used [25] . 

A cut of X = (14.33 GeV^/m^) = 0.59 has been suggested by the CLEO collaboration in 
order to eliminate the resonance region [8] . For this value the partially integrated rate is 

B„,, = 3.8 + 1.9 (4 + ^) - 134.7 + + 



2>m\ 



+614.9^ + 134.7^ + 560.2 A- (33) 
ml ml ml 

At this value of the cut x the coefficients of the nonperturbative matrix elements clearly in- 
dicate a poorly converging OPE. One can estimate the uncertainty induced by the 0(l/mf) 
parameters by fixing Aj to the values given in (27), then randomly varying the magnitudes 
of the parameters Pi.,% and /i between — (0.5GeV)^ and (0.5 GcV)^ as suggested by di- 
mensional analysis. We also impose positivity of pi as indicated by the vacuum saturation 
approximation [22] , and we enforce the constraint 



f^)'^ ^° Mli^MB{Mn + A) - MlAMpj MB + A) 



P2-%-% = ^^^ "^—r—zjm:-^^—^ (34) 



derived from the ground state meson mass splittings AMh = M^* — Mh {H = B,D) 
[18]. Here f3o is the usual coefficient of the beta function po = 11 — 2/3n/- = 25/3 for 4 
light flavours. Taking the 1 a deviation as a reasonable estimate of the uncertainties from 
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0{l/m\) contributions, we find the uncertainty in jBo.59 to be at the 10% level. It is clear 
from (33) that the pi contribution is large, and relaxing the positivity constraint on pi 
enlarges the uncertainty to about 20%. A similar statement can be made regarding the D0 
analysis [8] where the phase space cut is slightly higher, and the nonperturbative corrections 
are correspondingly somewhat larger. Since the cut on s cannot be lowered because of the 
■0' resonance, these uncertainties are intrinsic to our approach in the large dilepton invariant 
mass region. 

It is important to notice that in the invariant mass region below the J/ip resonance, 
the uncertainties from these matrix elements are much smaller. For example, integrating 
the differential spectrum up to the cut specified in the CLEO analysis [5] s — {Mjj^ — 
O.lGeVy/ml = 0.35 we find 

/ ds— = 22.0 1 + 0.5 4 + ' , +1-2 A + 7 , -3.74 - 1-2-^ ■ 
Jo.oi ds ml m% m^ 6m% ) m% m% ) 

(35) 

It is still true that the coefficient of the p\ term is ~ 10 times larger than that of the Ai 
term, but both sets of nonperturbative corrections are small relative to the parton level 
result in this region. Although this does not allow us to draw a strong conclusion about the 
convergence of the OPE, we can conclude that in this region the 0(l/m^) nonperturbative 
corrections are not a significant source of theoretical uncertainty. 



IV. THE FORWARD-BACKWARD ASYMMETRY 

The differential forward-backward asymmetry is defined by 

dA /•! , dB /-o , dB 



ds 

where 



[' dz^- C dz (36) 
Jo dzds J -I dzds 



z — cos9 — , ^ (37) 

u{s,ms) 

parameterizes the angle between the b quark and the in the dilepton CM frame. It has 
been shown [4] that new physics can modify this spectrum, so it is interesting to see how 
0{l/ml) terms contribute to the SM prediction. 

Integrating the triple differential decay rate (16) we find 

dA_ I «n 4(3 + 2g + 3g^)Ai , 4 (7 + lOg - 9g^) A^ 
lE-^' C,o 1^-8(1-.) + 

^4_(5 + 2s + s^) pi 4 (7 + 10s - Qs^) p2 \ 



3ml ^6 



^r^^Sf^^r ( 2s (3 + 2s + 3s^) A^ , 2s (9 + 14s - 15s^) A^ 

+ C, (s)6io l^-4s(l - s) + 

2s(l + 2s + 5s^)pi _ 2s(l + 6s-15s^)p2 ^ ..on 

3mf ml ) ^"^^^ 
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Here we have again omitted the trivial dependence on 7i_4 . 

It is clear from this expression that the third order terms do not have abnormally large 
coefficients, and therefore introduce only small variations relative to the second order ex- 
pression. 

An experimentally more useful quantity is the normalized FB asymmetry defined by 

dA _ dA dB 

ds ds ds 

Unfortunately, this normalized asymmetry has inherited the poor behavior of the differen- 
tial branching ratio in the endpoint region. In Fig. 3 we illustrate the uncertainties of the 



dA 
d§ 

-0.5 

-1 
-1.5 

-2 

0.2 0.4 0.6 0.8 

t 

FIG. 3. The normalized forward backward asymmetry. The three curves show the mean value 
and the la uncertainty of the forward backward asymmetry, obtained in a way similar to that 
explained in section III. 

normalized FB asymmetry originating from the matrix elements of the dimension six opera- 
tors. The three curves show the mean value and the la uncertainty of the forward backward 
asymmetry, obtained in a way similar to that explained in section III. We can see that up to 
a value of s = 0.7 the uncertainties are small, but for larger values of the dilepton invariant 
mass the uncertainties increase rapidly. Because of the necessity of the cut to eliminate 
the cc resonances, the accessible high dilepton invariant mass region is therefore restricted 
to a few hundred MeV. For the dilepton invariant mass region below the J/il) resonance 
{s < 0.35) the uncertainties are small. 




V. EXTRACTING A AND Ai FROM THE HADRON INVARIANT MASS 

MOMENTS 

Throughout this paper we have fixed the values of the leading non-perturbative param- 
eters A, Ai, and A2. However, these values must be determined from experiment. The 
most sensitive observablcs for this purpose are those which vanish in the parton model. It 
is interesting to ask how severely our ignorance of the values of the 0{l/ml) parameters 
compromises our ability to extract the values of A and Ai from such a measurement. It 
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has been suggested by Ali and Hiller [10] that one use the first two moments of the hadron 
invariant mass spectrum defined by 

(^h)= j^S^-^nT^dSH. (40) 

This idea is similar to the approaches used in the semileptonic B Xciv^ [18] and the rare 
radiative B — > decays [17], though the experimental task is more difficult in this case 
due to the small size of the branching ratio. To calculate these hadronic moments we relate 
them to calculable partonic moments via 

„ A(Ai + 3A2) 
^^""^-^ Mb 

+ fM|-2MBA + A^ + A, + 3A2-^i±^ + ^^ + ^^ + ^^^^ + ^^ 



y 2Mb 2Mb 

Ti + Ta + 3(T2 





(41) 



2Mb 

{S]i) = (m% - AMlK + 6M|A2 + 2M|(Ai + 3A2) - 4MbA3 - 4MbA(Ai + 3A2) 

-Mb(pi + 3p2) + Mb{Ti + T3 + 3(T2 + %))) {sD 
+4 (m^A' - 2MbA3 - MbA(Ai + 3A2)) {xD 

+ (4M|A - 12M2 A^ - 2M|(Ai + 3A2) + 12MbA3 + 10MbA(Ai + 3A2) 

+Mb(pi + 3p2) - Mb{Ti + Ta + 3(r2 + %))) (xoSo) 
+2 (m| A2 - 2MbA3 - (Ai + 3 A2) Mb A) (sq) 

+4MbA3(xo) (42) 
where we have used the mass relation 

appropriate at this order in the OPE. We therefore have to calculate the first two moments 
of the parton energy {xo), (xq) and parton invariant mass (sq), (sq), as well as the mixed 
moment {xoSq). Defining 



M(-'") = = ^ r^'"'^ dxo r; dso <s--^, (44) 



we give the results for the required partonic moments in Appendix C. As before, we have 
included the dependence on the cut on the lepton invariant mass in these results. It is im- 
portant to note that the results for the partonic moments given in Appendix C are expressed 
in terms of the b quark mass rrib and must be re-expressed in terms of the B meson mass 
Mb using the mass relation (43). Using again for the cut on the invariant mass the value 
proposed by CLEO > 14.33 GeV^, we find for the two moments 
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{Sh) = Ml 



0.36^ + 0.64-^ + 0.67-^ - 0.09^ + 8.48^ + 3.79 



+1.09 



'Ml 



Ml 



Ml 



Ml 



AAi 



4.88 



.AA. 



0.4l£r + 0.73^i±^ + 0.31 



Ml 



M% 



Ml 

% + 3r4 

Ml 



(45) 



M| 



_0.05^ + 0.14 



M| 



Ml 



0.53 



Pi 
M| 



0.21^ + 0.63^ + 0.46-"^^^ 



Ml 



Ml 



Ml 



-0.05 



Ti + 3r2 

Ml 



(46) 



Consider first the expression for {8^). The Ai term has a small coefficient and tends to 
cancel against higher order corrections, making this moment particularly insensitive to Ai. 
We can see the problem another way by solving this equation for Ai: the solution exhibits 
a pole near A = 0.4, close to the expected value of A [12]. As a result, the extracted value 
of Ai is extremely sensitive to the values of the higher order parameters. Since the presence 
of this pole persists as the value of the cut is changed, we conclude that this observable is 
unsuitable for extracting Ai. 

For the first moment {Sh) the convergence of the OPE is much better. Estimating the 
uncertainties from the unknown values of the dimension six operators by the method ex- 
plained in section III, we present the resulting constraint in the A — Ai plane in Fig. 4. 
Superimposed on this figure is the ellipse obtained in [27] from an analogous study of mo- 



-0.1 
-0.2 
-0.3 
-0.4 
-0.5 



0.2 0.4 ^ 0.6 0.8 

FIG. 4. The constraint in the A — Ai plane from {Sh)- The width of the band is entirely due 

to uncertainties from ©(l/m^) operators. The elHpse is the equivalent constraint from moments 
of semileptonic B Xcii^i- It is only the relative orientation, and not location, of the constraints 
which has meaning. 

ments oi B Xcii'i. Unfortunately, the bound from our analysis is nearly parallel to the 
major axis of this ellipse and, since it is only the relative orientation of the constraints which 
has meaning in this figure, this moment does not provide much additional information about 
the values of A or Ai. 
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VI. CONCLUSIONS 



Our purpose in this paper has been to study nonperturbative uncertainties in the rare 
inclusive decay B — > Xgi'^i'. Building on previous studies which evaluated the leading 
nonperturbative corrections, we have parameterized the corrections arising at 0{l/ml) in 
terms of two matrix elements of local operators pi 2, four matrix elements of non local 
operators Tl_4, and one matrix element of a four fcrmion operator fi. 

The numerical values of these parameters are unknown, yet even so a knowledge of 
the analytic form of the corrections allows us to study the convergence properties of the 
operator product expansion in various regions of phase space. Furthermore, the assumption 
that these parameters, being of nonperturbative origin, should be 0{AqQY)) permits us to 
make numerical estimates of theoretical uncertainties in observable quantities. 

We first considered the corrections to the differential spectrum dB/ds. The experimental 
spectrum contains two prominent resonances due to intermediate J/ip and ip' production, 
and the necessity of cutting these resonances out divides the accessible spectrum into two 
parts: the region of low dilepton invariant mass below the J/ip resonance, and the region 
of high dilepton invariant mass above the ip' resonance. In the first region, we find that 
the parton level calculation dominates, and that nonperturbative corrections are small. The 
operator product expansion appears to be converging according to expectation and, should it 
be possible to take experimental data in this region, the results will not suffer from significant 
nonperturbative uncertainties. 

On the other hand we do find that, as expected, the nonperturbative uncertainties in- 
crease as one moves into the high dilepton invariant mass region. It is well known that as 
one approaches this endpoint, the operator product expansion breaks down. The interesting 
result of our study, however, is that the expansion breaks down somewhat earlier than was 
anticipated, with the Oll/mf) uncertainties coming to dominate the integrated rate once 
the available range of s is reduced to about one quarter of the full range. We also showed 
that the rate obtained by integrating over the entire region above the ip' resonance contains 
uncertainties from dimension six matrix elements at the 10% level, the uncertainties being 
dominated by those from the pi matrix element. This result may impact the potential for 
doing precise searches for new physics using data from this region of phase space. 

We also studied the contributions from dimension six operators to the forward-backward 
asymmetry. This quantity probes different combinations of Wilson coefficients than the 
rate, and has been proposed as a compUmcntary source of information about possible new 
physics effects. As with the rate, the spectrum dA/ds is divided into two regions by the 
J/'^ and 'ip' resonances. We find that the dimension six contributions are not unduly large 
anywhere in the phase space, suggesting that this observable has a well behaved OPE. 
However if the differential forward-backward asymmetry is normalized to the differential rate, 
the resulting spectrum contains large uncertainties in the high dilepton invariant mass region. 
In the region of phase space below the J/%1) resonance, however, we find the nonperturbative 
corrections to be small. 

Finally we addressed a recent proposal suggesting that hadron invariant mass moments 
of the differential spectrum for B Xs£^£^ could, due to the sensitivity of these moments 
to nonperturbative effects, be used to constrain the values of the HQET parameters Ai and 
A. In the low Sh region the moments are suppressed relative to the rate and, considering the 
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already tiny branching ratio, it is unlikely that experimental measurements in this region 
will be forthcoming. Therefore we focused our attention on the high invariant mass region. 
In this region, we found that the first of these moments (Sh) provided a constraint in the 
A — Ai plane, but that this constraint was nearly the same as those derived from other, more 
experimentally promising, processes, and therefore seems to be of limited interest for this 
purpose. As for the second invariant mass moment (Sfj), we found that the nonperturbative 
uncertainties were such that it was not possible to extract a stable constraint on the values 
of A or Ai. From these results, we conclude that these moments are not well suited to the 
extraction of these parameters. 
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APPENDIX A: THE CONTRIBUTION OF DIMENSION SIX OPERATORS TO 

THE FORM FACTORS 

In this appendix we present the form factors Tj. These form factors have been calculated 
previously up to 0{l/ml) [9], and we do not reproduce those results here. We decompose 
the new contributions arising at C(l/m^) as 

For completeness we have included the full rhg dependence in these expressions, though in 
our analysis we set = 0. Defining x = 1 — 2v ■ q + s + ie with q = —, s = <f and 
Pi = Pi/iTT-h ; *® t^^^ third order contributions are 

- ^ ((2 + s) pi - 3 (2 - s) p2 + (pi + 3P2) v-q-{p, + 3p2) {v ■ qf) 
~3x^ (Pi + 3P2) il-v-q)(^s-{v- qf) 

+ ^Ji{l-v-qf{s-{vqf) (A2) 

Tf^^''"^' = --^(/5i + 3p2) 
bx 



(4pi + 6p2 - {pi + 3p2) V ■ q) 



3x2 
2 

+ (3sp2 - 2 (2pi + 3p2) v-q + 2 (2pi + 8^2) {v ■ qf) 
+±-^p,(l-vq){s-{vqf) (A3) 
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2 

+— ^ (l-v-q) (3p2 - (pi + 3p2) • q) 

T^f ' = ((Pi + 3p2) ((l - 5m^) s-2{l + ml){v qf)) 

-^T^ (s ((4m2 + s + m^s _ 4) + 3 (^8m^ + s + m^s) ps) 

+ (1 - 5m^) s (pi + 3p2) ^; • g + (1 + ml) ((8 - s) pi - Ssps) • qf 
-2 (l + m^) (Pi + 3p2) (vqf) 

~3i^ (1 + m^) (1 - • q) {s (pi + 3p2) - 2 (2pi + Sps) v ■ q) (s - {v qf) 
(s + 3mp + (1 + ^s) s ■ g - 2 (1 + ml) {v ■ qf) (A5) 



T^'^:^{l + ml) (pi + 3p2) 



3sx 

3s,T 



^T7^ (1 + ) (4P1 - (Pi + 3P2) ■ q) 



(1 + m^) (3sp2 + 2 (2pi + 3p2) v-q-2 (2pi + 3p2) • qf) 



- (1 + ^il) Pi{l-v-q)(s-{v qf) (A6) 

^^^^' = -^(1-^^) (Pi + 3p2)^-g 
2 

- 0--2 2 {'^-'rhl)v-q ((8 - s) pi - 3sp2 - 2 (pi + 3p2) v ■ q) 

(s (2pi + 3p2) + HPi + 3P2) ■ g - 2 (2pi + 3p2) {v ■ qf) 
(1 -ml){l-v q) pi(s^ -2svq- s{v ■ qf + 2{v ■ qf) (A7) 



3s2a;4 



2) 

X 



--^ (2 (3 + m^) s (pi + 3p2) - 6 (4p2 + s (pi + 3p2)) • q 

(j S X 

+2 (3 - mt) (pi + 3p2) {v ■ qf) 
■-^ ((1 - rhfj pi - 3m^p2) {1 - v q) (s - {v ■ qf) 

(j S X 

^ - pi {1 — V ■ q) [s^ — (1 — m^) sv ■ q — s{v ■ qf 



3^ A I 
SX^ 

+ (1 - m^) . qf) (A8) 
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yC.(C.±c,o) = ^ (p, + 3p2) + —,P2 {fnl + v q) (A9) 

*J X Ob 

(j S X 

(^1 + ^'^sPi + Sm^ps) (1 - v • g) V • g 
+77^^ (l + ^9 Pi (1 - ^ • g) - • g)') (AlO) 

(j S X 

APPENDIX B: THE DILEPTON INVARIANT MASS SPECTRUM WITH FULL 

MASS DEPENDENCE 



In this Appendix we present the dilepton invariant mass spectrum for a finite s-quark 
mass nis- The spectrum originating from operators of dimension d <b has been presented in 
Eq. (47) of [9]. The contributions from the time-ordered operators '7i_4 can be obtained by 
making the replacement (13) in this equation. Since the dilepton invariant mass distribution 
is independent of the definition of the four velocity of the heavy quark, the contribution 
proportional to p2 is related by reparameterization invariance [28] to the A2 contribution. It 
can be obtained by the replacement 



P2_ 

mi,' 



in the results of [9] . 

Thus, the only term we have to add to the existing literature to obtain the complete 
expression including all l/m\ contributions is the term originating from the Darwin operator 
whose matrix element is pi. This contribution is given by 



dB. 



pi 



ds 



Bopi 



licl 



(1^1 



v?{s,rhs) 



+ h5{su - s) 



(Bl) 



(B2) 



(B3) 



The function 'u(s, m^) = y (s — (1 — 'ms)'^){s — (1 + rhs)"^) is singular at the upper limit 
of integration. To regulate this divergence we defined a "star function" 

[F(s, rhs)], = hm | F(s, m,)^(s„ - s - (5) - 5{su - s - (3) r~ds F{s, m,) I . (B4) 

/3^0 I Jst I 



with limits of integration defined in (17) 

si< s < Su, 

with 

si^4rhf, Su^{l-rhs)^. 
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This "star function" is analogous to the common plus distribution. 
The functions appearing above are 



[1 - sf (ll + 13 s - 9 + 23 + 10 s 

- (50 + 37s + 48s^ + 38s^ + 70s^ + 45 s^) 
+ (85 + 150 s + 216 + 178 + 75 s^) - 2 (30 + 69 s + 72 + 25 s 



3^ ml 



+ (5 + 19 s) 7h° + (14 + 15 s) rh;^ -5 7% 



.12 



(B5) 



7^V^ = - 



8 

9i 



[l-sf (22-7s + 9s'^ + 19s^ + 5s^) 



- (78 - 34 s - 105 + 376 - 40 + 18 - 5 s 

+ f 70 + 67 s - 222 - 188 - 32 - 15 s^) mt 



6^ 



jCrCg _ 



+ (50 + 60s + 258 5^ + 136 s^) ml - (llO + 157s + 111 - 505^) ml 
+ (38 + 2s-75s2) + 9 (2 + 5s) -lOm]^ 

l-s)^(3-s + 17s^-3s=^) 
(lO - 13 s + 56 + 58 - 10 - 5 s^) 



(B6) 



+ (5 - 22 s + 12 - 34 - 25 s^) + 2 (lO + 29 s + 36 f + 25 s^) ml 

- (35 + 67 s + 50 s^) ml + (22 + 25 s) - 5 (B7) 

/, = -15^^(c?.+ (cfW+2Cf)') (B8) 
3^1 -4m2 V ^ ^ J 

Notice that these terms correctly reproduce the expression (26) in the limit rhg 0. 

APPENDIX C: THE MOMENTS UP TO ©(l/Ml) WITH A CUT ON THE 

DILEPTON INVARIANT MASS 

We write the moments in the form 

^(m,n) _ |o fcloMi^^{l> + Icfl^Mi'J}'''^ + M<S;^'"^ + cfMi'^'^'A , (CI) 
V ' ' ' ' / 

rr 

where is given in Eq. (32). The coefficient Cg depends on the parameters Xq and Sq as 
explained in section II, so we express the moments Mq^'""^ and M^*^'"^ as integrals which we 
evaluate numerically, 
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(m,n) 



9,9 — 



1\ rla-x) 



16 

,2 

2 



^Xoj-SOy^(C9±Cio)2 



■\cfixo,So)\' 



- So lm\ 



(m,n) 
7,9 



16 / 1\ rii^-x) f4 



For the other contributions we find 

\4 



^/xl^^Im< 



-'"10,10 — 



M. 



(1,0) 
7,7 



M- 



(2,0) 
10,10 



M 



(2,0) 
7,7 



30 



mt 



\2 (2{l-xfx (l + 6x-4x')^ 



mt 



pi /67- 30x + 30x'-35x'-32x^ 



45 



^ P2 / l + 10x-50x' + 55x'-16x' ^ 



mi 



2 (41 - 60 X + 18 x' + 4 - 3 x' + 24 log(x)) 

9 

/ 8(8-9x + x' + 61og(x)) \ 
ml\ 9 ) 

A2 /4 (7-2x-10x' + 5x' + 12 1og(x))\ 



+— 
mi 



+ 



pi /8 (2-27x + 19x' + 6x' + 18 1og(x))' 



27 



^P2_ ( 8 (13- 27x + 23x' -9x^-6 log(x))^ 



9 



2(1 - 2xo + so)^^'^^"^' 



2(l-2xo + so)rf'^^'^^"'^ 

(C2) 



^ (1-x)' (4 + 5x) ^ / (1-x)' (43 + 67x + 25x^) \ 



45 



270 



^ / (1-x)' (13 + 24 X - 222 x' + 125 x') ' 



mt 



90 



pi /(1-x)' (24- 27x- 3x' + 191x' + ir5x')' 



mt 



270 



P2 /(1-x)' (14- 63x- 306x^ + 175 x^')^ 



90 



2 (119 - 210 X + 120 x^ - 20 x=^ - 15x" + 6 x' + 60 log(x)) 

45 

Ai / 127 - 150 X - 12 x' + 44 x' - 3 x' - 6 x' + 84 log(x) ^ 



mt 



27 



(C3) 



(C4) 



(C5) 
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A2 / 9-2x + 4x^-36x=^ + 35x'-10x' + 121og(x) \ 
ml[ 3 ; 

Pj_ / 127 - 210 X - 24 x^ + 188 - 39 x' - 42 + 60 log(x) \ 
ml [ 27 ; 

^P^ / 39 - 66x - 8x^ + 84 x^ - 63 x' + 14x' + 12 log(x) \ .^^g. 

^io,i) _ Ai / (1-xf (13 + 19x + 8x^) \ 
' [ 30 ) 

^ X, gi-xf (3+13x-8x^) \ 
ml\ 6 / 

pi / (1-x)' (177 + 254x + 201x^ + 88x') \ 
ml\ 90 J 

^P2 ( {l-xf (l-18x-47x^ + 24x^) ^^ 

^ I 10 ) ^ ^ 

Ai / 2 (23-12x-18x^ + 4x' + 3x' + 24 1og(x)) \ 

^6 V 9 y 

A2 ( 2 (31-28x-18x^ + 20x=^-5x' + 24 1og(x)) A 

3 ) 

Pi / 2 (77 - 132 X - 5 1 x^ + "6 + 33 x ' 120 log(x)) \ 
rnl \ 27 J 

^ p (281 - 324 X - 54 x^ + 124 x^ - 27 x' + 168 log(x)) \ 
ml\ 9 / 



(0,1) 



M- 



(0,2) _ 

10,10 — ~~~2 



Ai /^4(l-x)' (4 + 5x)^ 



pi / 2(l-xr (31 + 64x + 40x^) \ 
ml \ 135 y 

+4 (2(l-x)' (l-26x + 10x-)\ 

mf \ 45 y ^ ^ 

(0,2) ^ A^ / 8 (119 - 210 X + 120 x^ - 20 x' - 15 x' + 6 x' + 60 log(x)) \ 
^'^ V 135 ) 

pi / 8 (139 - 60x- 210x^ + 140x' + 15x^-24x^ + 120 log(x)) \ 
ml \ 135 / 

P2 /16 (73 - 165x + 165x^-100x^ + 30x^-3x^ + 30 1og(x))\ 
ml\ 45 y 

_ Ai / (l-x)' (23 + 62x + 50x^) \ 
' [ 270 ) 

/ (1-x)' (13 + 82x- 50x^) \ 
ml V 90 J 



(CIO) 
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Pi Al-x)' (71 + 183 x + 276x^ + 190 X^) \ 
ml \ 270 j 

_P2_ Al - x)' (13 + 9 X - 252 x^ + 110 x') \ 
ml\ 90 ) 

M) ^ _ A / 4 (2 + 15x-33x^ + 16x=^ + 3x'-3x' + 6 1og(x)) \ 
rnl\ 27 j 

/ 4 (10 - 13x - x' + 8x' - 5x' + X' + 6 \og{x)) \ 
ml\ 3 ; 

^ pi 4 (62 - 75 X - 195 x' + 280 x' - 15x' - 57 x' - 30 log(x)) 



(Cll) 



ml \ 135 

P2 (A (257 - 285 X - 195 x^ + 400 x^ - 210 x^ + 33 x^ + 150 log(x)) ^ 

3 
b 



mM 45 ^ ' 
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